T 




Europaisches Patentamt 
European Patent Office 
Office europ6en des brevets 



<Ti) Publication number : 0 588 594 A3 



EUROPEAN PATENT APPLICATION 



<2i) Application number : 93307225.8 
Date of filing : 14.09.93 



<§) int CI. 5 : G05B 13/02, B01D 3/42 



(So) Priority : 14.09.92 US 944645 



© Date of publication of application : 
23.03.94 Bulletin 94/12 

(84) Designated Contracting States : 
BE DE FR GB IT 

@ Date of deferred publication of search report : 
24.08.94 Bulletin 94/34 

<g) Applicant : TEXACO, Inc. 
3901 Briarpark 

Houston, Texas 77042-5301 (US) 

@ Applicant : N E U RAL WARE , Inc. 
Penn Center West IV-227 
Pittsburgh, Pennsylvania 15276 (US) 



Inventor : Graettinger, Timothy J. 

221 Bluestone Drive 

Bethel Park, Pennsylvania 15102 (US) 

Inventor : Dubose, Paul A. 

2013 S. Lakeshore Drive 

Chapel Hill, North Carolina 27514 (US) 

Inventor : Federowicz, Alexander J. 

1814 Shady Avenue 

Pittsburgh, Pennsylvania 15217 (US) 

Inventor : Bhat, Naveen V. 

2100 S. Gessner, Apt 1506 

Houston, Texas 77063 (US) 

Inventor : Braden, William B. 

2229 West Alabama 

Houston, Texas 77098 (US) 

Inventor : Heckendoorn, Kent E. 

2718 Cason 

Houston, Texas 77005 (US) 



ffi) Representative : B runner, Michael John 
GILL JENNINGS & EVERY 
Broadgate House 
7 Eldon Street 
London EC2M 7LH (GB) 



CO 

< 



in 

to 
oo 
m 



Q. 

UJ 



@ Control system using an adaptive neural network for target and path optimization for a multivariate, 
nonlinear process. 

variables. The process is repeated until an opti- 
mal path is obtained, at which time the manipu- 
lated variable values are applied to the actual 
process. On a periodic basis all of the distur- 
bance, manipulated and controlled variables 
are sampled to find areas where the training of 
the neural network is sparse or where high 
dynamic conditions are indicated. These values 
are added to the set of values used to train the 
neural network. 



57) A control system having four major com- 
ponents : a target optimizer, a path optimizer, a 
neural network adaptation controller and a 
neural network. In the target optimizer, the 
controlled variables are optimized to provide 
the most economically desirable outputs, sub- 
ject to operating constraints. Various manipu- 
lated variable and disturbance values are 
provided for modeling purposes. The neural 
network receives as inputs a plurality of settings 
for each manipulated and disturbance variable. 
For target optimization all the neural network 
input values are set equal to produce a steady 
state controlled variable value. The entire pro- 
cess is repeated with differing manipulated vari- 
able values until an optimal solution develops. 
The resulting target controlled and manipulated 
variable values are provided to the path 
optimizer to allow the manipulated variables to 
be adjusted to obtain the target output. Various 
manipulated variable values are developed to 
model moves from current to desired values. In 
this case trend indicating values of the manipu- 
lated and disturbance variables are provided to 
produce time varying values of the controlled 
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(g) Control system using an adaptive neural network for target and path optimization for a multivariable, 
nonlinear process. 



(57) A control system having four major com- 
ponents : a target optimizer, a path optimizer, a 
neural network adaptation controller and a 
neural network. In the target optimizer, the 
controlled variables are optimized to provide 
the most economically desirable outputs, sub- 
ject to operating constraints. Various manipu- 
lated variable and disturbance values are 
provided for modeling purposes. The neural 
network receives as inputs a plurality of settings 
for each manipulated and disturbance variable. 
For target optimization all the neural network 
input values are set equal to produce a steady 
state controlled variable value. The entire pro- 
cess is repeated with differing manipulated vari- 
able values until an optimal solution develops. 
The resulting target controlled and manipulated 
variable values are provided to the path 
optimizer to allow the manipulated variables to 
be adjusted to obtain the target output. Various 
manipulated variable values are developed to 
model moves from current to desired values. In 
this case trend indicating values of the manipu- 
lated and disturbance variables are provided to 
produce time varying values of the controlled 
variables. The process is repeated until an opti- 
mal path is obtained, at which time the manipu- 
lated variable values are applied to the actual 
process. On a periodic basis all of the disturb- 
ance, manipulated and controlled variables are 
sampled to find areas where the training of the 



neural network is sparse or where high dynamic 
conditions are indicated. These values are 
added to the set of values used to train the 
neural network. 
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The invention relates to control systems used in 
process control environments, particularly to control 
systems developing optimal target outputs and taking 
optimal paths from current to desired target condi- 
tions for non-linear processes, and more particularly 5 
to control systems utilizing neural networks in the op- 
timization processes. 

In many industrial environments there are numer- 
ous processes which must be controlled. Examples 
include oil refineries, chemical plants, power plants, io 
and so on. For example, in oil refineries base crude 
oil feed stocks enter the refineries and after numer- 
ous individual distillation or reaction processes, a 
panoply of hydrocarbon streams are output. Each one 
of t he individual processes must be basically control- 15 
led or disaster could result. Above this basic standard, 
each individual process must be properly controlled to 
develop acceptable outputs. 

Each of the processes has numerous inputs or 
disturbance variables, such as input compositions, 20 
feed rates and feed stock temperatures. Numerous 
outputs or controlled variables also usually are pres- 
ent, such as the properties of various output streams. 
To perform process control, certain settable items or 
manipulated variables are also present, such as heat- 25 
ing, cooling and recirculation rates. 

Conventionally, control of the process was per- 
formed using feedback techniques. The controlled va- 
riables or outputs were measured and compared 
against desired values. If errors were present, 30 
changes were made to the manipulated variables 
based on the errors. One problem generally resulting 
from feedback control was stability concerns in the 
process, as the time delay may be critical in physical 
systems. Alternatively, the output may have reflected 35 
an undesirable trend or condition too late for adequate 
control of the process. Historically, proportional-inte- 
gral-derivative (PID) controllers were utilized for feed- 
back control. One problem generally resulting from 
the use of PID controllers was lack of multivariate 40 
control, as a PID controller only has one input variable 
and one output variable. Further, no economic control 
can be readily performed. Most real world systems 
have numerous controlled, disturbance and manipu- 
lated variables. Very complex arrangements of PID 45 
controllers could be utilized, but the sheer complexity 
often limited confidence and testability, much less a 
determination of optimal control. 

Certain feedforward techniques were utilized 
with some success. These techniques utilized a linear so 
model of the process being controlled and various 
changes to the manipulated variables were simulat- 
ed. Preferably the simulation allowed a quicker re- 
sponse time than feedback control and the model pre- 
ferably utilized all the major disturbance, manipulated 55 
and controlled variables of the process, allowing mul- 
tiple variable control. By using the disturbance or in- 
put variables, the controlled variables could be ob- 
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tained, those generally being present downstream of 
where the disturbance variables were measured, 
thus leading to the feedforward definition. However, 
the existing feedforward, modeling techniques have 
a major limitation. They operate primarily on linear or 
linearized systems, where a linear model can be de- 
veloped. The problem is that most of the more com- 
plex processes are by nature non-linear. Further, the 
models are not updated when minor process changes 
occur. Thus accuracy and/or range are sacrificed 
when linear modeling is utilized. 

Therefore a control system which can perform 
fast feedforward control for non-linear systems with 
high accuracy over a wide range is desirable. 

The present invention is a control system which 
provides accurate feedforward control for non-linear 
systems. Optimization techniques are utilized both to 
determine the optimal target manipulated and con- 
trolled variable values based on the disturbance vari- 
able values, various variable limits and economic con- 
cerns and to determine the optimal path from existing 
manipulated variable settings to the necessary set- 
,. tings to obtain the target controlled variable values. A 
neural network is used to model the process being 
controlled, allowing non-linear operation. Further, the 
neural network is designed to provide either steady 
state outputs, for use in the target optimization, or 
next time step outputs, for use in the path optimiza- 
tion. 

A control system according to the present inven- 
tion has four major components: a target optimizer, a 
path optimizer, a neural network adaptation controller 
and a neural network utilized by all three. The two op- 
timizers utilize sequential quadratic programming 
(SQP) to determine the various values provided to 
the neural network for modeling purposes, preferably 
a variation which uses only feasible values in the va- 
rious iterations to a solution, referred to as FSQP. 

In the target optimizer, the controlled variables 
are optimized to provide the most economically desir- 
able outputs, subject to the various system con- 
straints. The sum of the output stream values less the 
sum of the operating costs is maximized. In other in- 
stances other properties could be maximized, de- 
pending on the particular system. The FSQP portion 
provides possible manipulated variable values. These 
values are utilized along with the disturbance variable 
values in the neural network to produce an output. 
The neural network is designed to receive as inputs 
a plurality of settings for each manipulated variable 
and each disturbance variable. In this case all the val- 
ues for each manipulated and disturbance variable 
are set equal to produce a steady state controlled va- 
riable value. The output values of the controlled vari- 
ables and the related costs associated with the ma- 
nipulated variables are processed to develop an eco- 
nomic value. The entire process is then repeated with 
new manipulated variable values provided by the 
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FSQP portion. This iterative process continues until 
an optimal solution develops. If none develops, cer- 
tain restraints as to the feasible limits are removed or 
changed and the process repeated, until the most op- 
timal solution is developed. 

The target controlled variable values and associ- 
ated manipulated variable values which provide the 
optimized result are provided to the path optimizer to 
allow the manipulated variables to be adjusted from 
their present values to the values needed to obtain 
the optimized output Again an optimal output is de- 
sired. In this case the sum of the squared errors be- 
tween the actual controlled variable value and the de- 
sired value for a series of time steps and a manipulat- 
ed variable move aggressiveness factor is minimized 
as the objective function. The FSQP portion receives 
the various manipulated variable limits and develops 
various manipulated variable values to move from 
current to desired values. These values are combined 
with the disturbance values in the neural network. In 
this case time varying values of the manipulated and 
disturbance values are provided to the plurality of ma- 
nipulated and disturbance variable inputs of the neu- 
ral network. This results in the output of the neural 
network being the value of the controlled variables at 
the next time step. The process is repeated for each 
set of manipulated variable values developed by the 
FSQP portion so that hopefully the target values are 
obtained. When reached, the objective function value 
is stored, along with the manipulated and controlled 
variable values for that particular path, to allow deter- 
mination of the gradient of the objective function 
based on the manipulated variable changes. Then 
the entire process is again repeated for another path 
between the current and desired values. After the op- 
timal path is obtained, the first manipulated variable 
values of the optimal sequence are applied to the ac- 
tual process. 

Preferably the path optimization process is re- 
peated on a periodic basis, but the target optimization 
process occurs only when certain trigger limits for the 
various values are exceeded. Thus target optimiza- 
tion is performed less frequently, reducing the 
changes to the path optimizer, resulting in a more sta- 
ble output and requiring significantly less computa- 
tional time. 

In addition, a system according to the present in- 
vention also adapts to changing conditions to fine 
tune its operation. On a periodic basis all of the dis- 
turbance, manipulated and controlled variables are 
sampled. If the sampled data indicates areas where 
the training of the neural network is sparse or high dy- 
namic conditions are indicated, the values are added 
to the set of values used to train the neural network. 
Eventually over time, a very complete training set is 
gathered to allow best development of the neural net- 
work coefficients. Preferably retraining of the neural 
network is accomplished on operator command, but 



the retraining could be initiated automatically upon 
ceipt of a certain number of new samples. Certain 
checks are performed on the new coefficients and 
then, after approval, the new neural network coeff i- 
5 cients are utilized and the process more accurately 
controlled based on the increased training data. 

The preferred embodiment utilizes a control sys- 
tem as summarized above to control a debutanizer in 
an oil refinery. Details are provided below. 
10 A better understanding of the present invention 

can be obtained when the following detailed descrip- 
tion of the preferred embodiment is considered in 
conjunction with the following drawings, in which: 
Figure 1 is a piping drawing of a depropanizer- 
15 debutanizer system controlled by a control sys- 

tem according to the present invention; 
Figure 2 is a block diagram showing the various 
inputs and outputs of a control system according 
to the present invention for use on the debutan- 
20 izerof Figure 1; 

Figure 3 is a block diagram illustrating the inputs 
and outputs of the target optimizer in the control 
system of Figure 2; 

Figure 4 is a block diagram illustrating the inputs 
25 and outputs of the path optimizer in the control 

system of Figure 2; 

Figure 5 is a block diagram illustrating the inputs 
and outputs of the adaptation control in the con- 
trol system of Figure 2; 

30 Figure 6 is a diagram of the neural network used 

in the target optimizer, path optimizer and adap- 
tation control of Figs. 3, 4 and 5; 
Figure 7 is a block diagram of the electrical circui- 
try according to the present invention for use on 

35 the debutanizer of Figure 1; 

Figure 8 is a flowchart of the operation of the tar- 
get optimizer of Figure 3; 

Figure 9 is a flowchart of the operation of the op- 
timization process in the target optimizer of Fig- 

40 ure 8; 

Figure 10 is a flowchart of the operation of the 
path optimizer of Figure 4; 
Figure 11 is a flowchart of the optimization proc- 
ess in the path optimizer of Figure 10; 

45 Figure 12 is a flowchart of operation of the adap- 

tation control of Figure 5; 

Figure 13 is a flowchart of operation of the re- 
training task of Figure 12; and 
Figure 14 is a flowchart of operation of a schedu- 
50 ler for starting the target and path optimizers and 

adaptation control of Figs. 8, 10 and 12. 
Referring now to Figure 1, the piping drawing of 
the physical system of the preferred embodiment is 
shown. The preferred embodiment is the process 
55 control of a conventional debutanizer 20 in an oil re- 
finery. A depropanizer 22 has a bottoms outlet pipe 
24. A temperature sensor 26, preferably a thermocou- 
ple, and a pressure transducer 28 measure the tern- 
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perature and pressure of the stream in the pipe 24. 
Additionally, an analyzer 30 measures the C3 and C4 
hydrocarbon components in the stream in the pipe. 
C3 components are considered as those hydrocarbon 
compositions in the stream having 3 carbon atoms, 5 
such as propane propalyne. C4 components are 
those compositions in the stream having 4 carbon 
atoms, such as butane, butalyne and butadiene. Sim- 
ilarly C5 components are those compositions in the 
stream having 5 carbon atoms, such as pentane, pen- 10 
talyne and pentadiene. Conventionally the analyzer 
30 is a gas chromatograph, but in the preferred em- 
bodiment it is preferably an analyzer based on an in- 
ference technique as described in Patent Serial No. 

, entitled "Method for Inferring Con- is 
centrations of Process System Components for use in 
Control Systems", filed August 25, 1992. In the infer- 
encing technique the material, which is near its bub- 
ble point in the depropanizer 22, has its temperature 
and pressure measured at one point and is transmit- 20 
ted across an expansion valve to allow the stream to 
partially flash and then has its temperature and pres- 
sure remeasured. Because part of the stream will 
have vaporized, so that the pressure and temperature 
will have changed, this measuring of two temperature 25 
and pressure sets provides sufficient data points so 
that a linear regression technique can be used to de- 
termine two of the three component percentages, 
with the third being developed by simple subtraction. 
This inferencing is preferred because it is significantly 30 
faster then a gas chromatograph. 

Further downstream on the pipe 24 is a flow me- 
ter 32. The pipe 24 then proceeds to a valve 34, the 
valve opening being controlled by a level controller 
based on the level of the fluids in the depropanizer 35 
22. A feed pipe 36 is provided from the valve 34 to the 
debutanizer 20, preferably entering the debutanizer 
20 halfway up the debutanizer 20. A temperature sen- 
sor 38, again preferably a thermocouple, and a pres- 
sure sensor 40 determine the temperature and pres- 40 
sure of the feed stream to the debutanizer 20. Prefer- 
ably the debutanizer 20 has a series of trays as com- 
monly utilized in such devices. The debutanizer sys- 
tem basically produces two output streams, one of 
which is the C4 product and one of which is the C5 45 
product. 

The C5 product is obtained from the bottom of the 
debutanizer 20. A bottoms pipe 42 is connected from 
the bottom of the debutanizer 20 to a bottom valve 44. 
The opening of the bottom valve 44 is controlled by 50 
the level of fluid in the debutanizer 20 by a level con- 
troller so that the fluid level does not get too low or too 
high. A pipe 46 is connected from the valve 44 to the 
next location in the refinery process which will utilize 
and operate on the C5 product. The stream provided 55 
in this pipe 46 is considered to be an output of the 
process. An analyzer 48 is connected to the pipe 46 
to determine the concentration of C4 material in the 
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C5 product stream. 

Additionally, certain material from the bottom of 
the debutanizer 20 is reheated or rebofled and fed 
back into the bottom of the debutanizer 20. A reboil 
exit pipe 50 exits near the bottom of the debutanizer 
20 and proceeds through a heat exchanger 52 to a re- 
boil return pipe 54, which returns to the debutanizer 
20 at a level above the reboil exit pipe 50. A line 56 
containing 50 pound steam is provided through a flow 
sensor 58 to a valve 60. The opening of the valve 60 
is controlled by a flow controller to a desired set point 
based on the flow rate as measured by the flow meter 
58. The output of the valve 60 is provided through the 
heat exchanger 58 to provide the heat used to reboil 
the material obtained from the debutanizer 20. 

The C4 and C3 products in the feed stream are 
lighter than the C5 product and are removed from the 
top of the debutanizer 20. A top pipe 62 exits the top 
of the debutanizer 20 and is connected to a heat ex- 
changer 64. Preferably an analyzer 66 which meas- 
ures the C3 and C5 product concentrations is con- 
nected to the pipe 62. Cooling water is transmitted 
through the heat exchanger 64 to condense *and re- 
duce the temperature of the materials exiting the top 
of the debutanizer 20. A pipe 67 is connected from the 
heat exchanger 64 to deliver the cooled stream to a 
surge tank 68. A vent pipe 70 is provided from the top 
of surge tank 68 to transmit the C3 components which 
are still vaporized after the sub-cooling in the heat ex- 
changer 64. A pipe 72 is connected from the bottom 
of the surge tank 68 to a pump 74, whose output is 
connected to a pipe 76. A temperature sensor 71 
measures the temperature of the stream in pipe 72. 

The pipe 76 tees, with one direction proceeding 
through a flow meter 78 to a reflux valve 80. The re- 
flux valve 80 is controlled by a flow controller to a de- 
sired set point based on the flow rate provided by the 
flow meter 78. A pipe 82 is connected to the valve 80 
and to near the top of the debutanizer 20 to allow re- 
fluxing of the stream as desired. The second portion 
of the tee of pipe 76 goes through a flow meter 84 to 
an exit valve 86. The exit valve 86 is a pressure con- 
trolled valve based on a pressure reading taken at ap- 
proximately the middle of the debutanizer 20. The 
output of the valve 86 is connected to a pipe 88, which 
then goes to the next item in the chain which receives 
the output C4 product material. This is a simplified ex- 
planation of the piping of a debutanizer. 

The various disturbance variables in the debutan- 
izer control stream include the feed stream flow rate 
as measured by the flow meter 32, the feed stream 
temperature as measured by the temperature sensor 
38, the C3 and C4 component feed concentrations as 
determined by the analyzer 30 and the temperature 
of the reflux stream as measured by the temperature 
sensor 71 connected to the pipe 72. The controlled 
variables of the process are the concentration of C5 
materials exiting the top of the debutanizer 20 as 
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measured by the analyzer 66 and the C4 concentra- 
tion in the C5 product stream provided in the pipe 46 
as determined by the analyzer 48. The manipulated 
variables are the reflux flow rate as determined by 
the flow meter 78 and set by the valve 80, the reboll 5 
steam flow rate as measured by the flow meter 58 
and set by the valve 60 and the pressure in the deb- 
utanizer tower 20 as controlled by the valve 86. It is 
noted that each of these valves 60, 80 and 86 has an 
associated controller, with -the controller set points 10 
being adjustable to perform the control operations of 
the process. 

Turning now to Figure 2, a basic block diagram of 
the control system for use with the debutanizer 20 of 
Figure 1 is shown. Control system 100 has three ma- 15 
jor portions, a target optimizer 102, a path optimizer 
104 and adaptation control 106. Inputs to the control 
system 1 00 are the disturbance variables, the control- 
led variables and a series of other parameters utilized 
in the optimization processes. Some of these miscel- 20 
laneous inputs include the values of certain quantities 
of the C4 and C5 products, preferably done as a cost 
curve or table; the cost of the 50 pound steam utilized 
in the reboiling operation; and the cost of the cooling 
water utilized in the heat exchanger 64. To control the 25 
optimization operations it is also understood that the 
various manipulated variables have certain physical 
limits, so reflux flow rate limits, reboil steam rate limits 
and tower pressure limits are provided as inputs to the 
control system 100. Further, for reasons explained 30 
below relating to the maximum step size to allow a 
smooth transition taken in a change, rate of change 
limits for the manipulated variables are inputs. Addi- 
tionally, there are also certain target limits for the out- 
put stream impurities, that is the C4 concentration in 35 
the bottom or C5 product output and the C5 concen- 
tration in the overhead or C4 product output, for the 
outputs to be considered optimal. Additionally, other 
inputs include disturbance variable trigger values as 
will be explained below relating to the various rates of 40 
change in the disturbance variables necessary to ini- 
tiate the target optimization process. 

The primary outputs of the control system 100 
are the three manipulated variable values, that is the 
reflux flow rate set point, the reboil steam flow rate set 45 
point and the tower pressure set point as utilized with 
the valves 60, 80 and 86. Additionally, the predicted 
manipulated variable values and trending of the con- 
trolled variables are providing to the system to allow 
an operator to view the predicted future operation of so 
the process. 

Addressing the target optimizer 102, Figure 3 
shows a block diagram illustrating the target optimizer 
102 with its various inputs and outputs. In addition, a 
neural network 1 08 is shown as being incorporated in 55 
the target optimizer 102. The target optimizer 102 re- 
ceives the actual controlled variable values as meas- 
ured by the analyzers 48 and 66 and the disturbance 
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variable values. Additionally, certain miscellaneous 
inputs including the various cost numbers as provided 
to the control system 100 are provided to the target 
optimizer 102, as well as the bottom C4 and overhead 
C5 concentration limits and the reflux flow rate, reboil 
flow rate and tower pressure limits. Additionally, rate 
of change limits forthe manipulated variables are also 
provided. This is done as the path optimizer 104 can 
move only a finite number of steps in its prediction, so 
it is desirable to limit the change in the manipulated 
variables as requested by the target optimizer 102 to 
a value equal to the rate of charge times the number 
of steps. The outputs of the target optimizer 102 are 
the target or desired manipulated variable values and 
the target or desired controlled variable values to pro- 
duce the optimal operation of the debutanizer 20. The 
development of these values will be explained in de- 
tail below. 

Figure 4 illustrates a block diagram of the path 
optimizer 104 with its various inputs and outputs. The 
path optimizer 104 receives as inputs the target ma- 
nipulated variable and controlled variable values from 
the target optimizer 102, the actual values of tjje dis- 
turbance variables and the actual values of the con- 
trolled variables. Also, the reflux flow rate, reboil 
stream flow rate and tower pressure limits are provid- 
ed as miscellaneous inputs. The path optimizer 104 
also utilizes or contains the neural network 108. The 
outputs of the path optimizer 104 are the actual ma- 
nipulated variable values for the next time step to be 
applied to the system of Figure 1 and the future val- 
ues for use in the graphical trend display provided on 
the operator workstation. 

The adaptation control 106 (Fig. 5) is utilized to 
determine the initial coefficients or values used in the 
neural network 108 and to further train the neural net- 
work 1 08 as operations proceed. The adaptation con- 
trol 106 receives the values of the disturbance vari- 
ables, the manipulated variables and the control va- 
riables on a periodic basis and provides as outputs 
new weighting coefficients or the weight matrix for 
use by the neural network 108. 

The neural network 108 is shown in detail in Fig- 
ure 6. It is preferably a conventional three layer net- 
work with an input layer, a hidden layer and an output 
layer. The neural network 108 is preferably developed 
using matrix mathematical techniques commonly 
used in programmed neural networks. Input vectors 
are multiplied by a weight matrix for each layer, the 
values in weight matrix representing the weightings or 
coefficients of the particular input to the result being 
provided by the related neuron. An output vector re- 
sults. The input layer uses the inputs to the neural 
network 108 as the input vector, and produces an out- 
put vector used as the input vector for the hidden lay- 
er. The hidden layer matrix multiplication in turn pro- 
duces an output vector used as the input vector for 
the output layer multiplication. The output vector of 
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the output layer is the final output of the neural net- 
work, in this case the C4 and C5 impurity concentra- 
tion values or controlled variables. Preferably the 
neurons in the neural network 108 use a hyperbolic 
transfer function such as 

e* - e-* 

e* + e~* 

for X values of -1 to 1. The neural network 108, based 
on its training, readily handles nonlinear situations, so 
that a model of a nonlinear process is readily devel- 
oped by the use of the neural network 108. 

When used with the debutanizer 20, the neural 
network 108 has 32 input neurons, 10 hidden neurons 
and 2 output neurons. The inputs to the neural net- 
work 108 are the 4 sets of the manipulated variables 
and 4 sets of the disturbance variables. Four sets of 
the manipulated variables are utilized so that either 
steady state output values can be developed for use 
with the target optimizer 102 or time sampled output 
values, i.e. the next step, can be developed for the 
path optimizer 104. Sets of the disturbance variables 
are used because in the preferred embodiment the 
process is also dependent on changes in the distur- 
bance variables, as well as the manipulated vari- 
ables. So sets are used to better model the trend or 
time dependent nature of the disturbance variables. 
The neural network 108 is considered as being fully 
connected in that each of the neurons of the hidden 
layer receives the outputs of all the input layer neu- 
rons and each of the neurons of the output layer re- 
ceives the outputs of all of the hidden layers. The out- 
puts of the neural network 108 are the controlled va- 
riable values, the overhead C5 concentration and the 
bottom C4 concentration, which are preferably less 
than 2 mole %. 

Preferably the control system 100 is developed 
using general purpose computer systems, such as va- 
rious workstation computers or others, and the ap- 
propriate sensors and output drivers as shown in Fig. 
7. The various measured values are provided from the 
various analyzers 30, 48 and 66 and sensors 32, 38 
and 71 over a communications network from the 
measurement site to a data gathering or instrumen- 
tation and control computer 90 as appropriate for the 
specific plant environment. The input data is then 
transferred from the instrumentation computer 90 to 
the actual workstation computer 92 executing the 
steps of the control process. Similarly, the output val- 
ues are provided by the workstation 92 executing the 
control process steps to the instrumentation comput- 
er 90 and then out over the communications network 
to the particular valve controllers 60, 80 and 86. The 
various variable limits, cost information and trigger 
values are generally provided by the system operator 
directly to the workstation 92. Additionally, the neural 
network 108 is also preferably developed using the 
appropriate program on the workstation 92. 

As digital computers are preferably utilized in 
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executing the process, certain flowcharts are thus 
appropriate to indicate the operation of the actual 
process of the control system 100. Referring now to 
Figure 8, the target optimizer sequence 200 com- 
5 mences at step 202 where the computer reads vari- 
ous disturbance and controlled variable values. After 
reading is completed, control proceeds to step 204 to 
determine if any of the disturbance variable or con- 
trolled variables values changed by an amount great- 

10 er than the appropriate trigger value, which values 
have been previously set by the system operator us- 
ing certain utility functions. In conventional systems 
the target optimizer sequence 200 operates on a per- 
iodic basis, even if there has been little or no change 

15 in the actual input stream. This results in needless 
repetition of the target optimization procedure and 
has a possible result of variations in the target values, 
so that hunting or instability may result In the prefer- 
red embodiment the optimization process is prefer- 

20 ably performed only if certain levels of change have 
occurred, reducing the amount of time spent comput- 
ing the results and reducing the number of changes 
to the path optimizer or after a relatively long time if 
no changes sufficient to trigger a process have occur- 

25 red. If none of the disturbance or controlled variables 
have changed sufficiently to exceed a trigger level, 
control proceeds to step 205 to determine if the back- 
up or deadman time has passed since the optimiza- 
tion process was performed. This time is relatively 

30 long, preferably many sample periods, and is done 
just for added security as to the reasonable nature of 
the operation. If the period has not passed, control 
proceeds to step 206, where the sequence 200 is 
paused or goes to sleep until restarted or awakened 

35 by a scheduler 500 (Fig. 14). After being awakened, 
control proceeds to step 202 to reanalyze the distur- 
bance and controlled variable values. 

If the change in one of the disturbance or control- 
led variable values had exceeded the trigger level as 

40 determined in step 204 or the long backup time had 
elapsed as determined in step 205, control proceeds 
to step 208 where the optimization procedure is per- 
formed. This is preferably done by utilizing the neural 
network 108 for modeling the process, a feasible se- 

45 quential quadratic programming technique to develop 
iteration values and maximization of the net value of 
the output streams. This operation is performed in a 
separate sequence 250 illustrated in Figure 9 and ex- 
plained below. After the sequence 250 completes, 

50 control proceeds to step 21 0 to determine if indeed an 
optimal solution had been obtained. In certain cases 
this may not be the case as the FSQP program may 
have timed out because it has become stuck in a local 
minima or cannot determine a transition in a case re- 

55 quiring a large move. If an optimal solution has not 
been obtained, control proceeds to step 212, where 
the various limits and constraints applied to the optim- 
ization process are relaxed. This restraint relaxation 
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can be developed utilizing one of four techniques. In 
one case, all of the various limits or clamps on the ma- 
nipulated variables and controlled variables are re- 
leased, so that the best solution that can be obtained 
is utilized. Alternatively, the limits of the least valuable 5 
products in terms of cost of the output product can be 
relaxed successively. That is, the limits of least valua- 
ble first, the next least valuable second and so on can 
be relaxed. In a third technique, a predefined order- 
ing, other than least valuable, can be developed and 10 
used to relax the order of the various constraints or 
limits. 

In the preferred embodiment a fourth technique 
is used where certain of the constraints are consid- 
ered soft constraints and others considered hard con- is 
straints. Hard constraints are typically based on phys- 
ical limits, such as steam rates and steam change 
rates, and so should not be changed. Soft constraints, 
on the other hand, are typically those not based on 
physical concerns but on other concerns. For exam- 20 
pie, the controlled variable limits are usually imposed 
but are not physical limits. Thus these soft constraints 
are more suitable to adjustment. If optimal solutions 
are not developed within the specified limits of the 
soft constraints, then in the preferred embodiment, 25 
the transfer curve, or value curve in the case of the 
preferred embodiment used in the optimization proc- 
ess is slightly altered so that the hard boundaries are 
somewhat relaxed, so that a transition to zero value 
for a given concentration occurs, not as a discontinu- do 
ity or step, but as a continuous function. For example, 
if the C4 concentration limits are set at 1 and 2 mole 
percent, then a line is developed between the value 
curve at the two limits and a near zero point at 0.9 and 
2.1 mole percent, with a simple cubic spline smooth- 35 
ing the transitions at all points of discontinuity. This 
slightly expanded range allows the optimization proc- 
ess sufficiently more room to operate to allow an op- 
timal solution to be developed and yet not be too far 
from the desired limits. 40 

After the particular values are relaxed, control re- 
turns to step 208 to determine if an optimal solution 
can be developed with the new limit values. 

If an optimal solution had been developed, control 
proceeds from step 210 to step 211, where the ma- 45 
nipulated and controlled variable target values are 
provided to the path optimizer 104. Control then pro- 
ceeds to step 206 to wait until it is time for the next 
sample. 

As noted above at step 208, the target optimizer 50 
sequence 200 calls the optimizer sequence 250 to 
perform the actual optimization operations for the tar- 
get optimizer 102. The optimizer sequence 250 be- 
gins at step 252, where the various economic factors 
of the output streams and the manipulated variables 55 
are determined. These values have preferably been 
previously provided to the computer by the system 
operator. Control then proceeds to step 254 where 



the limits or constraints are determined. These limits 
include the controlled variable limits and the manipu- 
lated variable limits. Control then proceeds to step 
256, the first step of the FSQP program. The FSQP 
program is preferably FSQP 2.4 available from Pro- 
fessor Andre L. Tits at the University of Maryland in 
College Park, Maryland. The FSQP program is a ser- 
ies of routines which perform conventional sequential 
quadratic programming <SQP) techniques with the 
modification that the iterations used in reaching the 
minimized or maximized solution are performed using 
only feasible values. This feasible value limitation is 
particularly preferred for use with constraint limited 
physical systems, so that impossible states are never 
considered. Two basic techniques are used in th$ 
FSQP program in developing iterate values. The first 
is the use of an Armijo type arc search. The second 
is the use of a nonmonotonic search along a straight 
line. The FSQP program includes the ability to call 
specific functions or sequences to develop needed 
values. In the preferred embodiment then the FSQP 
internally calls the neural network 108 and a module 
that performs the optimization equation as steps in its 
overall process of determining an optimal solution. 
For this reason the FSQP program is shown in Fig. 9 
by the broken line encompassing several steps. The 
description of these internal operations is believed to 
better illustrate the invention. 

In step 256 the FSQP program selects the next 
manipulated variable values to be utilized in the op- 
timization process. The manipulated variable values 
selected by the FSQP program in step 256, along with 
the current disturbance variable values, are provided 
as inputs to the neural network 108 in step 258. As 
previously noted, the neural network 108 receives 
four values of each manipulated variable and each 
disturbance variable. This is to allow the output of the 
neural network 108 to be a time dependent output 
prediction. Many processes, such as the operation of 
the debutanizer 20, are time dependent, that is, they 
have a certain inertia or time lag. Responses are not 
instantaneous, so trend and history information is 
necessary to predict the output of the operation at a 
time in the future. When the neural network 108 is 
trained, the use of the various trend indicating ma- 
nipulated and disturbance variable values, along with 
the current controlled variable values, allows this in- 
ertia or time lag to be built into the neural network 
coefficients. Thus when trend indicating manipulated 
and disturbance variable values are provided, a time 
dependent controlled variable value is produced. 
However, if the various values of manipulated and dis- 
turbance variables are the same, that is the values 
provided for each manipulated and disturbance vari- 
able trend selection are identical, a steady state con- 
dition is simulated and the output of the neural net- 
work 108 is the steady state values of the controlled 
variables. Thus, by properly training the neural net- 
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work 108 and providing varying or equal values, time 
dependent or steady state outputs can be obtained 
from a single neural network. In this particular in- 
stance of optimization sequence 250, the four values 
for each particular manipulated and disturbance van- s 
able are all the same because a steady state solution 
is desired. With these inputs, the neural network 108 
then produces the model controlled variable values. 

A correction factor value is added to these con- 
trolled variable values in step 260. The correction fac- 10 
tor is based on the difference between the actual con- 
trolled variable values as last measured by the ana- 
lyzers 48 and 64 or as developed by linear regression 
from the last measured values and the modelled val- 
ues as determined by the neural network 108 for that 15 
same time sample. This provides a correction factor 
for inaccuracies in the model as embodied in the neu- 
ral network 108. The corrected values and the values 
of the manipulated variables and for the controlled va- 
riables for this iteration are then used in an optimiza- 20 
tion calculation in step 262, which utilizes the feed 
stream flow rate and the output concentration values 
to develop the economic value of the output streams, 
which are summed to produce an output stream val- 
ue. The costs of the various materials utilized, such 25 
as the steam and the cooling water, are subtracted 
from the output stream value to determine a net eco- 
nomic value for this particular group of manipulated 
variable values. It is understood that this particular 
optimization calculation is particular to the preferred 30 
embodiment example and other goals or optimization 
values could be selected. 

If the optimization process is not complete, as de- 
termined by the FSQP program and evaluated in step 
264, control returns to step 256, where the FSQP pro- 35 
gram generates another set of manipulated variable 
values. This cycle continues until all of the iterations 
of the FSQP program have been developed and pre- 
ferably the target values for the controlled variable 
values have been developed, or the FSQP program 40 
indicates certain problems, such as a timeout or in- 
ability to develop a feasible iterate. Control then pro- 
ceeds to step 266. In step 266 a determination is 
made as to whether an optimal value has been devel- 
oped. If the solution developed is non-optimal, then 45 
a flag is set in step 268 for use in step 210. If it is op- 
timal, control proceeds to step 270, where the appro- 
priate, optimal manipulated variable values and con- 
trol variable values associated with the maximum 
economic value are provided to a predetermined lo- so 
cation for use by the path optimizer 104. Control then 
proceeds from steps 268 or 270 to step 272, which is 
a return to the calling sequence, in this case the tar- 
get optimizer sequence 200. 

The path optimizer sequence 300 (Fig. 10) per- 55 
forms its complete sequence of operations on a per- 
iodic basis, preferably every control cycle. The path 
optimizer sequence 300 commences operation in 



step 302, where the particular target values are re- 
trieved. Control then proceeds to step 304, where the 
FSQP program, the neural network 108 and a mini- 
mization sequence are called to determine the opti- 
mal manipulated variable values for the future time 
steps. After this optimization is performed in se- 
quence 350, as shown in detail in Figure 11, control 
proceeds to step 306 where the manipulated variable 
values for the next time step are provided to the low 
level controllers. In the case of the preferred embodi- 
ment, valves 60, 80 and 86 have their set points 
changed to the new values. Additionally at this point 
the manipulated variable values for the entire predi- 
cated sequence can be provided to an operator work- 
station to allow a graphic presentation to show the opr 
erator the planned moves of the path optimizer 106. 
Control then proceeds to step 308, where the path 
optimizer sequence 300 goes to sleep. As noted, pre- 
ferably the path optimizer 300 performs on a repeti- 
tive basis, to allow transfer of the optimal controlled 
variable values to allow a path to be developed as de- 
sired. This timing function is performed by the 
scheduler 500. When the proper time period has 
elapsed, the path optimizer sequence 300 is awak- 
ened by the scheduler 500 and control transfers from 
step 308 to step 302, where the path optimizer se- 
quence 300 continues its operation to provide the 
next manipulated variable step. 

The path optimizer sequence 350 (Fig. 11) com- 
mences at step 352, where the physical or feasible 
manipulated variable limits, the rate of change limits 
and the current manipulated, disturbance and con- 
trolled variable values are retrieved. Control then pro- 
ceeds to step 354, where the target manipulated va- 
riable and controlled variable values are retrieved. 
Control then proceeds to step 355, where the correc- 
tion factor of the controlled variables is calculated. 
The analyzers 48 and 66 periodically provide actual 
controlled variable values. These values are then 
used to adjust or correct the predicted values pro- 
duced by the neural network 1 08. Because in the pre- 
ferred embodiment the control cycle time is less than 
the time for the analyzers 48 and 66 to produce read- 
ings, intermediate readings or values are developed 
using a linear regression on the analyzer data. This 
allows the correction factor to be calculated on each 
control cycle. 

After the calculation, control proceeds to step 
356, where manipulated variable values are selected. 
The FSQP program in this instance selects entire ser- 
ies of manipulated value moves to bring the process 
from the current values to the target values. The rate 
of change of the manipulated variable limits are util- 
ized in this selection so that each proposed manipu- 
lated variable step does not exceed certain limits, 
thus keeping the controlled process from being too 
upset or destablized by a very large change in a ma- 
nipulated variable in a relatively short time. Thus each 
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of the manipulated variable values provided in the set 
is within the change limits from the previous value. At 
this step 356 the trend indication values of the ma- 
nipulated and disturbance variable values are devel- 
oped for each time step. Four values are developed 
for each manipulated and disturbance variable using 
a filter equation. A simple exponential filter is used to 
develop the trend values. In the filter the current fil- 
tered values are determined (after the first time step) 
by the equations 

OVfiter. = ClDlW^ ♦ (1 - "^actual. 

MV 19Uttm = a MV^m + <1 - a) MV BCtuali0 
where a is between 0 and 1, with a lower a value in- 
dicating more weighting on current values and a low 
a value indicates more weighting on historical factors. 
Four different a or history factor values are used to 
show trend information. Because each a value pro- 
vides a different historical weighting, this allows sim- 
ple development of trend information. The four a val- 
ues may be the same for all the manipulated and dis- 
turbance variables, but preferably different a value 
sets can be used with each variable to better track the 
time dependent nature of each variable. The neural 
network 108 then receives the various values in step 
358 and produces an output of the controlled variable 
values at the next time step. In fhis case the neural 
network 108 receives trending manipulated and dis- 
turbance variable values, not equal values, so that a 
next time step prediction output is obtained from the 
neural network 108 for each set of input variables. 
This allows development of a time-based path from 
one state to the next The neural network 108 thus 
provides the controlled variable values for each step 
in the set, so that the predicted output is available for 
the entire set to develop the values for the entire path. 
The predicted controlled variable values are provided 
to step 360 for correction. Control then proceeds to 
step 362, where the optimization calculation is per- 
formed. In this case, optimization is slightly more 
complex than the equivalent operation in step 262. In 
this particular case the goal is to develop a minimized 
balance between errors and an aggressive or stable 
transition to a low error condition. To easily obtain this 
value, the difference between the target and correct- 
ed controlled variable values is squared and multi- 
plied by a first weighting factor. As a second compo- 
nent, the proposed change in the manipulated vari- 
able value is multiplied by a second weighting factor. 
The two values are then added to form the objective 
function value. This operation is performed for each 
step in the path, with each addition value further be- 
ing summed to develop a total factor. The first weight- 
ing factor controls the impact of errors in the minimi- 
zation process, while the second controls the aggres- 
siveness of operation, which relates to the stability or 
smoothness of operation. The weighting factors are 
set by the operator. With this summed value being 



minimized, this is considered to be the desired path 
from the current location to the desired goal. Control 
then proceeds to step 364, to determine if a mini- 
mized value has been obtained. If not, control pro- 

5 ceeds to step 365 to determine if an iteration count 
has been exceeded or calculation time exceeded. If 
not, control proceeds to step 356. This loop continues 
until the predicted result is optimal. This may require 
a number of iterations. 

w After some period, the controlled variable values 

should converge or the iteration or time limits are ex- 
ceeded. In those cases, control proceeds from step 
364 or step 365 to step 367, where the entire set of 
manipulated and controlled variable values fort he op- 
ts timal path is moved to the system. All of these values . 
are used by the operator's workstation to show plan- 
ned movements. Control then proceeds to step 374, 
where the manipulated variable values for the next 
particular step in the process associated with the op- 

20 timal path as determined by the minimum summed 
value are set up to be provided as an output of the op- 
timize sequence 350. Control then proceeds to step 
376, a return to the calling sequence, the path optim- 
• izer sequence 200. 

25 The adaptation control 106 executes the adap- 

tation sequence 400 (Fig. 12). The adaptation se- 
quence 400 is shown in Figure 13 and commences at 
step 402, where the sequence 400 obtains all of the 
current controlled variable, disturbance variable and 

30 manipulated variable values for possible storage as a 
record. The manipulated variable values are the fil- 
tered values as used in the path optimizer 104 and 
are stored for each sampling interval so that historical 
data is developed for the manipulated variables as 

35 needed to allow full training of the neural network 
108. Control then proceeds to step 404 to determine 
if the values of the disturbance variables and control 
variables indicate that the system is operating in a 
new region in which data has not been fully obtained. 

40 Preferably the number of disturbance variables and 
control variables sets up the dimensions of an array, 
in the preferred example of Figure 1 a sue dimensional 
array. Each one of these dimensions or variables is 
then further broken down into sub-intervals across 

45 the feasible range of the variables. This then devel- 
ops a multiplicity of cells, zones or regions relating to 
particular operating conditions of the debutanizer 20. 
In step 404 the computer performing the adaptation 
sequence 400 selects the particular disturbance val- 

so ue variable and control variable values of this partic- 
ular sample record and indexes to that particular re- 
gion as developed by the mapping technique. If this 
is a new operating region where there are fewer then 
a given number of samples, it is considered appropri- 

55 ate to save these for retraining of the neural network 
1 08. Control thus proceeds to step 428. If this was not 
a new operating region, control proceeds from step 
404 to step 408 to determine if there is a very high dy- 

9 
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namic change in the disturbance or controlled vari- 
ables. If so, t his sample may be of interest because it 
has been determined through experience that these 
high transient conditions are particularly desirable in 
training the neural network 108 and control proceeds 
to step 406. If there is no high dynamic situation, con- 
trol proceeds to step 410 to determine if an operator 
request to include this particular record has been in- 
dicated. If so, control proceeds to step 406. If not, con- 
trol proceeds to step 412 to determine if any other 
miscellaneous rules which may be developed have 
indicated that this particular sample is to be savecj. If 
so, control proceeds to step 406. If not, control pro- 
ceeds to step 414 to determine if the operator has 
provided a retraining command. A retraining com- 
mand indicates that it is an appropriate time to use the 
collected sample set to retrain the neural network 
1 08. Preferably retraining is commenced only on com- 
mand, but retraining could be automatically initiated 
after a long period of time, such as one month. 

If it is time to train the neural network 1 08, control 
proceeds to step 415 to determine if enough new 
samples are present to merit retraining. As retraining 
takes considerable time and resources and retraining 
with only a few new samples will result in only a very 
minor change, a predetermined number of samples is 
desired before retraining. If enough are present, con- 
trol proceeds to step 416, where the system determi- 
nes the training age limit, that is the oldest of the 
training samples to be utilized. For instance, if the de- 
propanizer 20 were to have received particular dam- 
age to one of the internal trays on a given day, it would 
not be appropriate to train the neural network 108 
with data from prior to that date, because that data 
would not be an accurate indication of the operation 
of the debutanizer 20 at the present time. Therefore 
by setting the training age limit, only data under the 
current existing physical conditions of the system 
would be utilized in retraining the neural network 108. 
If automatic retraining initiation is being used, a de- 
fault age limit could be selected. After step 416, con- 
trol proceeds to step 418, where a separate retraining 
task 450 (Fig. 13) is initiated. As noted, the process 
is preferably running on a multitasking environment, 
so this time intensive operation is separated and done 
in parallel with the real time operations. Control then 
proceeds to step 424. Control would also proceed 
from step 414 to step 424 if a training command had 
not been received and from step 415 to step 424 if too 
few new samples were present. In step 424 the adap- 
tation sequence 400 pauses or goes to sleep until 
next awakened by the scheduler 500 to sample the 
various variable values. 

In step 406 the computer determines the partic- 
ular number of samples in the region as defined by 
the controlled variables and disturbance variables 
values. Control proceeds to step 426 to determine if 
that particular region is full of samples. If not, control 



proceeds to step 428, which is also where control has 
proceed from step 404, where all of the particular va- 
riable values are added to the training set with a time 
stamp to allow both replacement of the sample set 
5 and utilization based on the training age limits as set 
forth in step 416. If the region is considered full, oon- 
trol proceeds from step 426 to step 430, where the 
oldest, non-sacred training model sample record is re- 
moved. There may be certain highly distinctive sam- 
10 pie records which are desired to be permanent por- 
tions of the training set. If so, these values can be in- 
dicated as being sacred records and thus non-remov- 
able. Control proceeds from step 430 to step 428 
where the new sample is added, replacing the re- 
ts moved record. After step 428 control proceeds to step 
424 to wait until the next time to sample the variables. 

The actual retraining task 450 is illustrated in Fig- 
ure 1 3. The task 450 commences at step 452 where 
the neural network 108 is retrained using convention- 
20 al techniques, preferably back propagation techni- 
ques, as well known to those skilled in the art Control 
then proceeds to step 454 where a steady state gain 
change analysis is performed. Step inputs ar$ provid- 
ed to the old and new weight matrices and th£ result- 
25 ing output changes compared. If the change in the 
output exceeds certain limits or has actually changed 
direction, then the new weight matrix is suspect, as 
t hese values s hould not change radically between re- 
training intervals. If suspect as determined in step 
do 456, control proceeds to step 464. If the steady state 
gains appear to be within the limits, control proceeds 
to step 458 where the new weight matrix is provided 
for operator approval. Operator approval before im- 
plementation is preferred to allow all interested par- 
35 ties to become comfortable with the proposed weight 
matrix and resultant process operation. If not ap- 
proved in step 460, control proceeds to step 464. If ap- 
proved, in step 462 the new weight matrix values re- 
place the prior values and operation of the control 
40 system 1 00 then proceeds with the new values. Con- 
trol then proceeds to step 464, where the retraining 
task 450 terminates. 

Operation of the scheduler 500 is shown in Fig. 
14. The scheduler 500 awakens the target optimizer 
45 sequence 200, the path optimizer sequence 300 and 
t he adaption control sequence 400 on a regular basis 
to allow those operations to be performed. Operation 
of the scheduler 500 commences at step 502 where 
the target optimizer sequence 200 is awakened or re- 
50 started. Control proceeds to step 504, where control 
loops until the target optimizer sequence 200 has put 
itself back to sleep. Control then proceeds to step 
506, where the pat h optimizer sequence 300 is awak- 
ened. Control proceeds to step 508, where control re- 
55 mains until the path optimizer sequence 300 returns 
to sleep. Control then proceeds to steps 510 and 512 
where similar actions are taken relating to the adap- 
tation control sequence 400. Control then proceeds to 
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step 514, to determine if the time for the next sample 
or pass through the sequences has elapsed. If not, 
control remains at step 514 until the sample time ar- 
rives. If so, control returns to step 502 and the proc- 
ess repeats. s 

Preferably, during normal operations the target 
optimizer 102, the path optimizer 104 and the adap- 
tation control 106 are all executing sequentially as 
shown in Fig. 14. But the target optimizer 1 02 and the 
path optimizer 104 cannot be allowed to control the 10 
debutanizer 20 until satisfactory coefficients are de- 
veloped for the neural network 1 08. Preferably this; is 
done by running just the adaptation sequence 400 on 
the debutanizer 20 as it operates under the control of 
its prior control system, which is to be replaced by the 15 
control system 100. This allows actual, site-specific 
data to be readily gathered to develop a close original 
weight matrix. As a next step, the control system 100 
can be run in control but no write shadow mode, with 
the target optimizer 102 and the path optimizer 104 20 
outputs not actually applied to the controllers. A com- 
parison of the two control systems can then be per- 
formed for an extended period to gain confidence in 
the control system 1 00. Eventually the outputs are ac- 
tivated and a changeover occurs. Then the adapta- 25 
tion control 106 can develop more accurate training 
information, which over time will tune the neural net- 
work weight matrix to the optimal values for the par- 
ticular debutanizer 20. This additional training can 
readily increase the effective nonlinearity of the mod- 30 
el to better approach the nonlinearities of the actual 
system, resulting in greatly improved performance 
over much broader ranges than a linear model at- 
tempting to control the same system. 

It is understood that the debutanizer 20 is the 35 
preferred embodiment and has been provided as an 
example. The techniques and processes according to 
the present invention can be utilized in many other 
process control environments, particularly multivari- 
ate and more particularly nonlinear, and no limita- 40 
tions are necessarily intended by the detailed descrip- 
tion of debutanizer operation. Further, it is understood 
that other neural network arrangements can be used, 
depending on the particular process and environ- 
ment. Additionally, the number of manipulated, distur- 45 
bance and controlled variables, optimization goals 
and variable limits can be changed to suit the partic- 
ular process of interest. 

Having described the invention above, various 
modifications of the techniques, procedures, material so 
and equipment will be apparent to those in the art. It 
is intended that all. such variations within the scope 
and spirit of the appended claims be embraced there- 
by. 



Claims 

1. A control system for a process having a plurality 

* 

of disturbance variables as process input condi- 
tions, a plurality of manipulated variables as proc- 
ess control conditions and a plurality of controlled 
variables as process output conditions, the con- 
trol system comprising: 

means for periodically determining the 
values of the disturbance variables; 

means for optimizing the values of the 
controlled variables based on the determined dis- 
turbance variable values and providing target ma- 
nipulated variable values associated with said 
optimized controlled variable values; 

means receiving said target manipulated 
variable values for optimizing the movement from 
the existing manipulated variable values to said 
target manipulated variable values and providing 
new manipulated variable values; and 

a neural network coupled to said control- 
led variable optimizing means and said move- 
ment optimizing means and receiving manipulat- 
ed variable values and disturbance variable val- 
ues and producing the resulting controlled vari- 
able values representative of the process being 
controlled. 

2. The control system of claim 1, wherein said neu- 
ral network has a plurality of inputs of each ma- 
nipulated variable, wherein when said plurality of 
inputs for each manipulated variable receive 
equal values, the resulting controlled variable 
values are steady state values and when said 
plurality of inputs for each manipulated variable 
receive different values, the resulting controlled 
variable values are future values. 

3. The control system of claim 2, wherein said con- 
trolled variable optimizing means and said move- 
ment optimizing means include means for devel- 
oping a plurality of manipulated variable values 
for provision to said neural network to produce a 
plurality of controlled variable values for the re- 
spective optimization, wherein said controlled va- 
riable optimizing means provides equal values to 
said plurality of inputs of each manipulated vari- 
able and wherein said movement optimizing 
means provides time dependent, varying values 
to said plurality of inputs of each manipulated va- 
riable. 

4. The control system of claim 1, wherein said neu- 
ral network has a plurality of inputs of each dis- 
turbance variable, wherein when said plurality of 
inputs for each disturbance variable receive 
equal values, the resulting controlled variable 
values are steady state values and when said 
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plurality of inputs for each disturbance variable 
receive different values, the resulting controlled 
variable values are future values. 

5. The control system of claim 4, wherein said neu- 5 
ral network has a plurality of inputs of each ma- 
nipulated variable, wherein when said plurality of 
inputs for each manipulated variable receive 
equal values, the resulting controlled variable 
values are steady state values and when said w 
plurality of inputs for each manipulated variable 
receive different values, the resulting controlled 
variable values are future values. 

6. The control system of claim 5, wherein said con- 15 
trolled variable optimizing means and said move- 
ment optimizing means include means for devel- 
oping a plurality of manipulated variable and dis- 
turbance variable values for provision to said neu- 
ral network to produce a plurality of controlled va- 20 
riable values for the respective optimization, 
wherein said controlled variable optimizing 
means provides equal values to said plurality of 
inputs of each manipulated and disturbance vari- 
able and wherein said movement optimizing 25 
means provides time dependent, varying values 

to said plurality of inputs of each manipulated and 
disturbance variable. 

7. The control system of claim 1, wherein said con- 30 
trolled variable optimizing means includes means 

for setting limits of the controlled and manipulat- 
ed variables to be utilized during optimizing and 
means for altering said controlled variable limits 
if an optimal solution cannot be developed using 35 
said limits. 

8. The control system of claim 7, wherein said 
means for altering said limits removes all limits. 

40 

9. The control system of claim 7, wherein said 
means for altering said limits alters said limits in 
order of least economic effect. 

10. The control system of claim 7, wherein t he optim- 45 
ization of the values of the controlled variables 
utilizes a transfer function between measured 
values and values to be optimized and wherein 
said means for altering said limits alters said 
transfer function to allow for a continuous function so 
outside of said limits. 

11. The control system of claim 10, wherein said con- 
tinuous function rapidly ramps to a zero optimiz- 
ing value outside of said limits. 55 

1 2. The control system of claim 1 , wherein said move- 
ment optimizing means operates on a regular, 

12 



periodic basis, and wherein said controlled vari- 
able optimizing means includes means for setting 
disturbance variable or controlled variable 
change limits and operates when the change in a 
disturbance or controlled variable exceeds its re* 
spective limit. 

1 3. The control system of claim 1 , further comprising: 

means for monitoring the disturbance va- 
riable values, the manipulated variable values 
and the controlled variable values and storing 
said monitored variable values for use in training 
said neural network. 

14. The control system of claim 13, wherein said 
monitoring and storing means determines the re- 
spective disturbance and controlled variable val- 
ues and stores said monitored variable values if 
the number of stored monitored variable values 
within a predetermined region of said distur- 
bance and controlled variable values is below a 
predetermined number. 

« 

15. The control system of claim 13, wherein said 
monitoring and storing means determines the dy- 
namics of said disturbance and controlled vari- 
able values and stores said monitored variable 
values if said dynamics exceed predetermined 
limits. 

16. A method of controlling a process having a plur- 
ality of disturbance variables as process input 
conditions, a plurality of manipulated variables as 
process control conditions and a plurality of con- 
trolled variables as process output conditions, 
the method comprising the steps of: 

periodically determining the values of the 
disturbance variables; 

optimizing the values of the controlled va- 
riables based on the determined disturbance va- 
riable values and providing target manipulated 
variable values associated with said optimized 
controlled variable values; and 

optimizing the movement from the existing 
manipulated variable values to said target ma- 
nipulated variable values and providing new ma- 
nipulated variable values; and 

said steps of controlled variable optimizing 
and movement optimizing utilizing a neural net- 
work receiving manipulated variable values and 
disturbance variable values and producing the re- 
sulting controlled variable values representative 
of the process being controlled. 

17. The method of claim 16, wherein said neural net- 
work has a plurality of inputs of each manipulated 
variable, wherein when said plurality of inputs for 
each manipulated variable receive equal values, 
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the resulting controlled variable values are stea- 
dy state values and when said plurality of inputs 
for each manipulated variable receive different 
values, the resulting controlled variable values 
are future values. 5 

18. The method of claim 17, wherein said steps of 
controlled variable optimizing and movement op- 
timizing each include developing a plurality of 
manipulated variable values for provision to said 10 
neural network to produce a plurality of controlled 
variable values for the respective optimization, 
wherein said step of controlled variable optimiz- 
ing provides equal values to said plurality of in- 
puts of each manipulated variable and wherein 15 
said step of movement optimizing provides time 
dependent, varying values to said plurality of in- 
puts of each manipulated variable. 

19. The method of claim 16, wherein said neural net- 20 
work has a plurality of inputs of each disturbance 
variable, wherein when said plurality of inputs for 
each disturbance variable receive equal values, 

the resulting controlled variable values are stea- 
dy state values and when said plurality of inputs 25 
for each disturbance variable receive different 

1 

values, the resulting controlled variable values 
are future values. 

20. The method of claim 19, wherein said neural net- 30 
work has a plurality of inputs of each manipulated 
variable, wherein when said plurality of inputs for 
each manipulated variable receive equal values, 

the resulting controlled variable values are stea- 
dy state values and when said plurality of inputs 35 
for each manipulated variable receive different 
values, the resulting controlled variable values 
are future values. 

1 

21. The method of claim 20, wherein said steps of 40 
controlled variable optimizing and movement op- 
timizing include developing a plurality of manipu- 
lated variable and disturbance variable values for 
provision to said neural network to produce a 
plurality of controlled variable values for the re- 45 
spective optimization, wherein said step of con- 
trolled variable optimizing provides equal values 

to said plurality of inputs of each manipulated and 
disturbance variable and wherein said step of 
movement optimizing provides time dependent, so 
varying values to said plurality of inputs of each 
manipulated and disturbance variable. 

22. The method of claim 16, wherein said step of con- 
trolled variable optimizing includes setting limits 55 
of the controlled and manipulated variables to be 
utilized during optimizing and altering said con- 
trolled variable limits if an optimal solution cannot 



be developed using said limits. 

23. The method of claim 22, wherein said altering of 
said limits removes all limits. 

24. The method of claim 22, wherein said altering 
said limits alters said limits in order of least eco- 
nomic effect 

25. The method of claim 22, wherein the optimization 
of the values of the controlled variables utilizes a 
transfer function between measured values and 
values to be optimized and wherein said step of 
altering said limits alters said transfer function to 
allow for a continuous function outside of said lirrtr 
its. 

26. The method of claim 26, wherein said continuous 
function rapidly ramps to a zero optimizing value 
outside of said limits. 

27. The method of claim 16, wherein said step of 
movement optimizing operates on a regular, per- 
iodic basis, and wherein said step of controlled 
variable optimizing includes setting disturbance 
or controlled variable change limits and operating 
when the change in a disturbance or controlled 
variable exceeds its respective limit 

28. The method of claim 16, further comprising the 
step of: 

monitoring the disturbance variable val- 
ues, the manipulated variable values and the 
controlled variable values and storing said moni- 
tored variable values for use in training said neu- 
ral network. 

29. The method of claim 28, wherein said step of 
monitoring and storing variable values determi- 
nes the respective disturbance and controlled va- 
riable values and stores said monitored variable 
values if the number of stored monitored variable 
values within a predetermined region of said dis- 
turbance and controlled variable values is below 
a predetermined number. 

30. The method of claim 28, wherein said step of 
monitoring and storing variable values determi- 
nes the dynamics of said disturbance and control- 
led variable values and stores said monitored va- 
riable values if said dynamics exceed predeter- 
mined limits. 

31. A system for performing a distillation process in 
an oil refinery having a plurality of disturbance 
variables as process input conditions, a plurality 
of manipulated variables as process control con- 
ditions and a plurality of controlled variables as 
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process output conditions, the system compris- 
ing: 

a distillation column having a means for re- 
ceiving an input feed stream and providing over- 
head and bottom output streams; 5 

means for cooling said overhead output 
stream; 

means for controllably returning a portion 
of said cooled overhead output stream to said dis- 
tillation column, said rate of controllably returning 10 
being a manipulated variable; 

means for controllably maintaining \he 
pressure of materials in said distillation column, 
said level of controllably maintaining the pressure 
being a manipulated variable; 15 

means for controllably heating a portion of 
material in said distillation column, said rate of 
controllably heating being a manipulated vari- 
able; 

means for periodically determining the 20 
values of the disturbance variables; 

means for optimizing the values of the 
controlled variables based on the determined dis- 
turbance variable values and providing target ma- 
nipulated variable values associated with said 25 
optimized controlled variable values; 

means receiving said target manipulated 
variable values for optimizing the movement from 
the existing manipulated variable values to said 
target manipulated variable values and providing do 
new manipulated variable values; and 

a neural network coupled to said control- 
led variable optimizing means and said move- 
ment optimizing means and receiving manipulat- 
ed variable values and disturbance variable val- 35 
ues and producing the resulting controlled vari- 
able values representative of the distillation proc- 
ess being controlled. 

32. The system of claim 31, wherein said neural net- 40 
work has a plurality of inputs of each manipulated 
variable, wherein when said plurality of inputs for 
each manipulated variable receive equal values, 

the resulting controlled variable values are stea- 
dy state values and when said plurality of inputs 45 
for each manipulated variable receive different 
values, the resulting controlled variable values 
are future values. 

33. The system of claim 32, wherein said controlled so 
variable optimizing means and said movement 
optimizing means include means for developing a 
plurality of manipulated variable values for provi- 
sion to said neural network to produce a plurality 

of controlled variable values for the respective 55 
optimization, wherein said controlled variable op- 
timizing means provides equal values to said 
plurality of inputs of each manipulated variable 



and wherein said movement optimizing means 
provides time dependent, varying values to said 
plurality of inputs of each manipulated variable. 

34. The control system of claim 31 , wherein said neu- 
ral network has a plurality of inputs of each dis- 
turbance variable, wherein when said plurality of 
inputs for each disturbance variable receive 
equal values, the resulting controlled variable 
values are steady state values and when said 
plurality of inputs for each disturbance variable 
receive different values, the resulting controlled 
variable values are future values. 

35. The control system of claim 34, wherein said neu- 
ral network has a plurality of inputs of each ma- 
nipulated variable, wherein when said plurality of 
inputs for each manipulated variable receive 
equal values, the resulting controlled variable 
values are steady state values and when said 
plurality of inputs for each manipulated variable 
receive different values, the resulting controlled 
variable values are future values. 

36. The control system of claim 35, wherein said con- 
trolled variable optimizing means and said move- 
ment optimizing means include means for devel- 
oping a plurality of manipulated variable and dis- 
turbance variable values for provision to said neu- 
ral network to produce a plurality of controlled va- 
riable values for the respective optimization, 
wherein said controlled variable optimizing 
means provides equal values to said plurality of 
inputs of each manipulated and disturbance vari- 
able and wherein said movement optimizing 
means provides time dependent, varying values 
to said plurality of inputs of each manipulated and 
disturbance variable. 

37. The system of claim 31, wherein said controlled 
variable optimizing means includes means for 
setting limits of the controlled and manipulated 
variables to be utilized during optimizing and 
means for altering said controlled variable limits 
if an optimal solution cannot be developed using 
said limits. 

38. The system of claim 37, wherein said means for 
altering said limits removes all limits. 

39. The system of claim 37, wherein said means for 
altering said limits are alters said limits in order of 
least economic effect. 

40. The control system of claim 37, wherein the op- 
timization of the values of the controlled variables 
utilizes a transfer function between measured 
values and values to be optimized and wherein 
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said means for altering said limits alters said 
transfer function to allow for a continuous function 
outside of said limits. 

41. The control system of claim 40 f wherein said con- s 
tinuous function rapidly ramps to a zero optimiz- 
ing value outside of said limits. 

42. The system of claim 31, wherein said movement 
optimizing means operates on a regular, periodic 10 
basis, and wherein said controlled variable optim- 
izing means includes means for setting distur- 
bance variable change limits and operates when 

the change in a disturbance variable exceeds its 
respective limit 1 s 

43. The system of claim 31, further comprising: 

means for monitoring the disturbance va- 
riable values, the manipulated variable values 
and the controlled variable values and storing 20 
said monitored variable values for use in training 
said neural network. 

44. The system of claim 43, wherein said monitoring 

and storing means determines the respective dis- 25 
turbance and controlled variable values and 
stores said monitored variable values if the num- 
ber of stored monitored variable values within a 
predetermined region of said disturbance and 
controlled variable values is below a predeter- 30 
mined number. 

45. The system of claim 43, wherein said monitoring 
and storing means determines the dynamics of 
said disturbance and controlled variable values 35 
and stores said monitored variable values if said 
dynamics exceed predetermined limits. 

46. The system of claim 31, wherein said means for 
controllably heating includes: 40 

means for receiving materials from ap- 
proximately the bottom of said distillation column; 
means for heating said received materials; 

and 

means for returning said heated received as 
materials to said distillation column. 

47. A method for developing a sample set for training 
a neural network, the method comprising the 
steps of: so 

obtaining values of the various inputs and 
outputs of the neural network at a specific time 
to form a new sample; 

developing an n-dimensional matrix of 
cells, wherein n is the total number of inputs and 65 
outputs of the neural network, each axis has a to- 
tal range corresponding to the value range of the 
respective input or output of the neural network, 



each axis total range being subdividing into cell 
ranges to result in a plurality of cell ranges for 
each axis, whereby the total number of cells in 
said matrix is the product of number of cell ranges 
in the total range for each axis; 

determining the corresponding cell in said 
n-dimensional matrix based on the obtained input 
and output values in said new sample; 

determining the number of previously stor- 
ed samples in said corresponding cell; and 

adding said new sample to the previously 
stored samples if said determined number is be- 
low a predetermined limit. 

48. The method of claim 47, further comprising the 
steps of: 

if said determined number of previously 
stored samples exceeds said predetermined 
number, determining if said new sample includes 
certain values indicating other predetermined 
conditions desirable for training the neural net- 
work; and 

removing a previously stored sarpple in 
said corresponding all if said new sample in- 
cludes said certain values and adding said new 
sample to the previously stored samples. 

49. The method of claim 48, wherein the previously 
stored sample which is removed is the oldest 
sample. 

50. The method of claim 47, wherein the previously 
stored sample which is removed is the oldest 
sample which is not marked as being a perma- 
nent sample. 



15 



BP 0 588 594 A2 




EP 0 588 594 A2 





3 3 



^ gf UJ 111 

g |8 8 

less 
c tea a 



If 



c s 

~1 



ee. 



CO 



CO 



D 
til 
© 



I 



CO 

cc 

UJ 

s 



Ul 

o 



o 



17 



EP 0 588 594 A2 




18 



€P 0 588 594 A2 






83 

So 
1" 



I 
I 

s 




ii 

> Ui 



O 
O 

I 



19 



EP 0 588 594 A2 




20 



EP 0 588 594 A2 



NEURAL NETWORK 



REFLUX FLOW RATE 



REBOIL STEAM FLOW RATE 



TOWER PRESSURE 



FEED FLOW RATE 



FEED TEMPERATURE 



FEED C3 CONCENTRATION 



FEEO C4 CONCENTRATION 



REFLUX TEMPERATURE 




106 



OVERHEAD 06 
CONCENTRATION 

BOTTOMS 04 
CONCENTRATION 



FIG. 6 



21 



EP 0 588 594 A2 



• 




CD 



22 



EP 0 588 594 A2 




FIG. 8 

23 



EP 0 588 594 A2 



{OPTIMIZE*) 280 

determine! 262 
economic 

FACTORS 



DETERMINE 
OUTPUT 
LIMITS 



2B4 



SELECT NEXT 
MVS 



256 



I 



EXECUTE 1 288 
NEURAL 
NETWORK 



I 



APPLY 
CORRECTION 
FACTOR 



i 



260 



CALCULATE I 282 
OPTIMIZATION 
EQUATION 





SET TARGET 
MVS A CVS 



270 



268 



1 



SET ERROR 
FLAG 



[ 



RETURN l 272 



PSQP 



FIG. 9 



24 



EP 0 588 594 A2 



PATH 1 300 
OPTIMIZER 



RECEIVE 
TARGET MVS 
& CVS 



302 




RUN OPTIMIZE 

o Determine 

MV6 



304 



PROVIDE NEXT I 306 
MVS TO 
SYSTEM 



308 



FIG. 10 

25 



EP 0 588 594 A2 




I DETERMINE MV 362 
AND MV CHANGE 
LIMITS AND DV, 
MVaCV VALUES! 



DETERMINE CV 
AND MV 
TARGET 




SET NEXT STEp] 

MVS FROM 
OPTIMAL PATH 



calculate] 395 

CORRECTION * 
FACTOR 



MOVE OVER MV 
ft CV SETS OF 
VALUES 



36? 



SELECT NEXT SET 
OF MV MOVES ft 
CALCULATE MV ft , «« 
DV TREND VALUES 1 356 




364 



REDICTED CVS > 
TARGET CVS? 



TERATION COUNT 01 
TIME EXCEEDED? 



FIG. 1 1 



26 



EP 0 588 594 A2 




27 



EP 0 588 594 A2 



RETRAIN 



I 



D 



450 



RETRAIN 
NEURAL 
NETWORK 



452 



I 



STEADY STATE 
GAIN ANALYSIS 



454 




PROVIDE 
MATRIX FOR 
APPROVAL 



456 




REPLACE 

WEIGHT 

MATRIX 



462 



C 




END 



464 



FIG 



28 



€P 0 588 594 A2 



r 




FIG. 14 



29 



